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Abstract. - We establish the qualitative behavior of the incommensurability e, optimal 
domain wall filling v and chemical potential /i for increasing doping by a systematic slave- 
boson study of an array of vertical stripes separated by up to d = 11 lattice constants. Our 
findings obtained in the Hubbard model with the next-nearest neighbor hopping t' = — O.lSt 
agree qualitatively with the experimental data for the cuprates in the doping regime x < 1/8. 
It is found that t' modifies the optimal filling v and triggers the crossover to the diagonal (1,1) 
spiral phase at increasing doping, stabilized already at a; ~ 0.09 for t' — — 0.3t. 



Interest in transition metal oxides has never been restricted to the most spectacular phe- 
nomenon of the high-Tc superconductivity, but also concerns, inter alia, metal- insulator tran- 
sitions, colossal magnetoresistance, and orbital ordering [1]. Among these phenomena the 
so-called stripe phases attract much attention. In this context the recent indications of uni- 
versal magnetic excitations in doped cuprates are intriguing [2]. Remarkably, some features 
of the magnetic spectra including their anisotropic two-dimensional character established in 
a detvi^inned YBa2Cu306+5 (YBCO) sample [3] can be understood in terms of fluctuating 
stripes suggesting that they are one of generic properties of the copper oxides [4]. 

Typically, such states result from the competition between the superexchange interaction, 
which stabilize the antiferromagnetic (AF) long-range order in the parent Mott insulator, 
and the kinetic energy of doped holes. Indeed, the magnetic energy is gained when electrons 
occupy the neighboring sites and their spins coupled by the superexchange order as in the 
Neel state, whereas the kinetic energy is gained when the holes can move and the AF order is 
locally suppressed along a domain wall (DW), leading to the formation of site-centered (SC) 
stripes. Other possible structures are bond-centered (BC) stripes. In this case, AF domains 
with a lower hole density and a stronger spin polarization are separated by DWs given by 
ladders with an increased hole density and a weak ferromagnetic order on the rungs. 

The most direct evidence for stripe phases in doped antiferromagnets has come from neu- 
tron scattering studies in which charge and spin modulations are identified by the appearance 
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of several extra incommensurate (IC) Bragg peaks [5]. Indeed, neutron diffraction measure- 
ments performed on Lai.6-rENdo.4Sr2;Cu04 (Nd-LSCO), a model compound for which the 
evidence of vertical spin and charge stripe order is strongest, revealed that magnetic peaks 
are displaced from the AF maximum at Qaf = (tTjTt) to the points Q = 7r(l ± 2e, 1) and 
Q = 7r(l, 1 ± 2e) [6]. Remarkably, the incommensurability e varies linearly with doping e = x 
in the underdoped regime of x < 1/8 meaning a fixed stripe fiUing v = 1/2. A similar value 
of the number of holes per DW v = 0.59 has also been estimated in a recent resonant soft 
X-ray scattering study of La2-a:Baa;Cu04 (LBCO) which is a more direct evidence of charge 
modulation [7]. In contrast, beyond x = 1/8, one finds in experiment a lock-in effect with 
e = 1/8, corresponding to a robust stripe phase with a charge (magnetic) unit cell consisting of 
four (eight) sites, and the AF domains with three atoms along the x direction. The essentially 
identical modulation and doping dependence of e was observed in superconducting crystals of 
La2-a:Sra:Cu04 (LSCO) with X > 0.05 [8]. Also in YBCO the incommensurability increases 
first linearly and next saturates at the effective hole doping x ~ 0.10 [9]. 

Apart from the neutron scattering, stripe phases have also measurable consequences in 
angle-resolved photoemission spectroscopy (ARPES) [10]. Indeed, the low-energy spectral 
weight of Nd-LSCO at a; = 1/8, is mostly concentrated along the antinodal T — X and T — Y 
directions, while there is a distinct gap for charge excitations around the S = (7r/2,7r/2) 
point, as expected in the SC stripe picture [11-13]. In contrast, ARPES spectra of both 
LSCO and Nd-LSCO at a; = 0.15 have revealed the existence of appreciable spectral weight 
along the nodal T — S direction. Therefore, as the BC stripes reproduce quite well the nodal 
segments [12], it seems that upon increasing doping they are formed at the expense of the 
SC ones. The relevance of the BC stripes at the doping level x = 0.15 is also supported by 
recent studies which have yielded pronounced spectral weight both in the nodal and antinodal 
directions, reproducing quite well the experimental results in Nd-LSCO and LSCO [14]. 

An alternative scenario which might explain the IC spin structure is a deformation of the 
AF order which optimizes a hole motion within a spiral phase [15]. The interest in the spiral 
state was recently renewed due to experimental results indicating a spin glass behavior of 
LSCO at small doping consistent with this scenario [16]. On the one hand, observation of 
the charge order in Nd-LSCO and LBCO rules out the spiral phase since the density of holes 
is then expected to be uniform. On the other hand, the question whether the charge order 
is a generic feature of the cuprates is far from being resolved yet [5]. In fact, Lindgard [17] 
has shown that spiral states can also resolve the universality of magnetic excitations in the 
cuprates and provide a competing paradigm with the stripe phase concept. 

Hence, an interplay between domain wall structures and spiral phases poses an interesting 
problem to study [18] even though the searching for the optimal filling of domain walls corre- 
sponding to the true ground state is rather a formidable task. Therefore, we limit ourselves 
to vertical site-centered (VSC) and vertical bond-centered (VBC) stripe phases, with the size 
of the AF domains varying from d = 11 to 3 lattice constants. In both structures, the largest 
distance d = 11 corresponds to a unit cell with 22 atoms. Note that the same length of the 
magnetic unit cell for a fixed d makes the SC and BC structures practically indistinguishable 
from each other in neutron diffraction experiments. 

We study the stripe and spiral phases using the Hubbard model. 



where the electron hopping t involves the nearest neighbor sites {ij), t' the next-nearest 
neighbor sites {{ij)), and U stands for the on-site Coulomb interaction. There are several 
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Fig. 1 - (Color online) Free energy gain 5F per site in the VSC stripe phases with respect to the AF 
phase as a function of doping x, as obtained for the t-t'-U model (1) with U — 12t and: (a) t' = 0; 
(b) t' = —0.15t. Domain walls are separated by d = 3, . . . , 11 lattice constants. Circles and squares 
show the corresponding data for (1,0) and (1,1) spiral order, respectively. 

Fig. 2 - (Color online) Free energy gain SF per site as in fig.Q obtained for the t-t'-U model (1) with 
U = 12t and t' — —0.3t for: (a) VSC stripe phases; (b) VBC stripe phases. 



experimental and theoretical studies suggesting that t' is finite in the cuprates. For example, 
topology of the Fermi surface seen by ARPES [10] can be only understood by introducing 
t' [19]. It also offers an explanation for the variation of Tc among different families of hole- 
doped cuprates [20]. Moreover, exact diagonalization (ED) studies have shown that while 
the d-wave superconductivity correlation is slightly suppressed by < in the underdoped 
region, it is substantially enhanced in the optimally doped and overdoped regions, indicating 
that t' is of great importance for the pairing instability [21]. 

We employ a rotationally invariant version of the slave-boson (SB) approach in spin 
space [22], which opens a possibility of studying both noncanted stripe phases and spiral 
order on equal footing. More details on the calculation method can be found in ref. [23]. Cal- 
culations were carried out on square clusters with the linear dimension along the x direction 
(perpendicular to the DWs) chosen as an even multiplicity of the elementary stripe unit cell 
dimension in all cases, so that the cluster size used for the considered unit cells varied from 
128x 128 up to 144x 144. These calculations became possible by developing an efficient scheme 
in the reciprocal space which makes use of the stripe symmetry [24]. For so large systems the 
finite size eflFects are below 10~^t and one is able to obtain a realistic comparison of the free 
energies of structures with different size of the unit cell. In the numerical studies, we have 
chosen U = 12t, which gives the ratio of J/t — 1/3 (with the superexchange J = At'^/U), 
being a value representative for LSCO [25]. 
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We begin with the free energy gain of the VSC stripe phase with respect to the AF 
phase, SF = F — Faf, shown in fig. ^as a function of doping x for representative values 
of the next-nearest neighbor hopping, i.e., t' = 0, t' = — O.lSt (VBC phases are addressed 
below) . Quite generally, one observes that the energy favors the stripe phases with the largest 
distance d = 11 between DWs for the lowest values of x, and the curves for next smaller d—1 
systematically cross the ones for d upon increasing doping, meaning that stripe phases with 
gradually decreasing AF domains become optimal phases. For t' — 0, this effect continues 
until the d — A stripe phase is reached as the d — 3 stripe phase is a highly-excited state in 
this case. On the contrary, finite t' results in a deeper energy minimum of the d = 3 phase 
and stabiHzes it in the overdoped regime x > 0.2. Next, except for the smallest d < 4 cases, 
increasing \t'\ shifts the free energy minima towards lower doping level which should affect 
the filling of the DWs and the charge distribution in the stripe ground state. Remarkably, for 
t' = — O.lSt, a value very close to that derived for LSCO [26], one finds that the most stable 
stripes are separated by d = 4 lattice spacings in a sizeable doping range above x ~ 1/8, both 
for the SC and BC stripe phase, in agreement with the neutron scattering experiment [8] and 
with dynamical mean field theory (DMFT) calculations [13] for LSCO. 

Regarding the spiral phases, upon hole doping away from half-filling the AF order becomes 
immediately unstable towards a diagonal (1,1) spiral order with Q = 7r(l — 2e, 1 — 2e) and 
then at higher doping x ~ 0.18 towards a vertical (1,0) one characterized by Q = 7r(l — 2e, 1), 
as suggested by early SB studies of the Hubbard model [27]. However, a finite next-nearest 
neighbor t' = —0.15t has severe consequences for the interplay between both spiral phases and 
shifts the crossover towards a higher doping x ~ 0.23 in agreement with the SB studies of the 
t-t'-J model [28]. It also affects the competition between the spiral and stripe phases so that 
at X > 1/8 the latter become unstable with respect to the (1,1) spirals. However, the energy 
difference between the most stable stripe phase with d = 4 and the lowest energy spiral phase 
is less than 0.002i, which is considered to be too small to lead to an unambiguous conclusion. 
Indeed, comparison of the SB data with available ED results on a 4 x 4 cluster obtained in 
the Hubbard model aX t' = and a; = 1/8 doping [29] shows that taking into account stripe 
phases allows one to approach closer the ED energy than the best spiral phase but due to the 
quantum fiuctuations the energy difference remains at the 9% level (c/. table Pi . 

In contrast, a larger value of the ratio \t' /t\ = 0.3 as expected for YBCO [26], clearly 
drives the system towards the diagonal (1,1) spiral phase and the crossover from the VSC 
(VBC) stripe phase appears already sHghtly below (above) x = 0.09, respectively, as depicted 
in fig. 13 Note that in the case of the SC DWs it is the d — A stripe phase which is unstable 
towards the spiral order and indeed the smallest distance between DWs established in YBCO 
corresponds to d = 5 [9]. Our results are then consistent with both density matrix renor- 
maHzation group [30] and Hartree-Fock [31] studies indicating the suppression of the stripe 
phases with increasing \t'\. Moreover, large removes also the degeneracy between the two 
stripe structures which shows up at t' = — 0.15t. Indeed, as reported in tabled the energy 
difference between the lowest energy SC and BC configuration in the doping regime x < 0.2 is 

Table I - Comparison of the free energy F per site for various phases as found within the SB 
approximation in the Hubbard model at x — 1/8 with the ED data of ref. [29]. The ground state (VBC 
and VSC) stripe phases are separated by d — 7 lattice constants. Parameters: U — 12t, t' = 0. 



phase 


AF 


(1,0) 


(1,1) 


VBC 


VSC 


ED [29] 


F/t 


-0.5393 


-0.5613 


-0.5700 


-0.5756 


-0.5756 


-0.6282 
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comparable to the accuracy of the present calculations. In contrast, when t' — — 0.3t, the SC 
stripe phase is noticeably more stable than its BC counterpart even for x < 0.2. The crossover 
towards the BC stripes above x = 0.2 remains in agreement with a recent analysis of the prop- 
agation of a hole inside the BC DW which has shown that such a DW structure enables a 
larger kinetic energy gain than a narrower SC one, which becomes especially important at a 
large doping level when the distances between stripes are small [32]. 

The experimental data should be compared with our findings concerning the doping 
dependence of the magnetic incommensurability e = l/2d and the optimal stripe filling, 
V = Nii/{NyN-£,w), where TVh is the hole excess compared to half-filling, whereas Ny stands 
for the actual length of the cluster with iVow DWs along the y direction. The ground state 
properties of both SC and BC phases are shown in fig.El The points in fig. |S|were deduced 
from figs, m and [21 and correspond to the middle of the stability region of the lowest energy 
configuration. The only exception is the d — 3 case with t' — —0.3t in which they are plotted 
for the minimum of the free energy. Such a choice guarantees that, at each particular doping 
level, a considered stripe phase with a given periodicity d would be indeed realized at least in 
the vast majority of the system. Comparing fig.|Slwith the experimental results in LSCO [7,8], 
one finds that the robust stability of the half-filled {i^ = 1/2) d — i stripes at x — 1/8 requires, 
in agreement with the previous SB studies on small 16x16 clusters [33], a finite next-nearest 
neighbor hopping t' < — 0.15i. In fact, for t' = —Q.l'ot present calculations give almost a 
linear dependence e = x for x < 1/8 and may be considered as reproducing a lock- in effect 
in a sizeable doping range above x ~ 1/8 until the d = 3 stripe phase sets in. We emphasize 
that stable d = i stripe phases with e = 1/6 have also been found in the SB studies of the 
three-band Hubbard model in the doping regime x > 0.225 [34]. The apparent absence of this 
phase in the experimental data for LSCO could follow from two effects which go beyond the 
present study: (z) repulsion between DWs, and (m) quantum fluctuations which destabilize 
the ladder- like AF domains. 

Remarkably, in the regime where e follows linearly x, an increasing density of stripes allows 

Table II - Comparison of the ground state free energy F per site for the VSC and VBC stripe phases 
as found in the t-t' -U model with U = 12t and: t' — — O.lSf and t' = — 0.3t. 







t' = 


-0.15t 






t' = 


~0.3t 








VSC 




VBC 




VSC 




VBC 


X 


d 


F/t 


d 


F/t 


d 


F/t 


d 


F/t 


0.050 


11 


-0.4263 


11 


-0.4263 


9 


-0.4289 


10 


-0.4280 


0.055 


10 


-0.4360 


10 


-0.4359 


8 


-0.4388 


9 


-0.4378 


0.060 


9 


-0.4456 


9 


-0.4455 


7 


-0.4488 


8 


-0.4477 


0.070 


8 


-0.4649 


8 


-0.4648 


6 


-0.4686 


7 


-0.4673 


0.080 


7 


-0.4841 


7 


-0.4840 


5 


-0.4882 


6 


-0.4869 


0.090 


6 


-0.5034 


6 


-0.5032 


5 


-0.5080 


5 


-0.5064 


0.100 


5 


-0.5225 


6 


-0.5224 


4 


-0.5275 


5 


-0.5260 


0.120 


5 


-0.5607 


5 


-0.5607 


4 


-0.5659 


4 


-0.5645 


0.140 


4 


-0.5985 


4 


-0.5983 


3 


-0.6022 


4 


-0.6005 


0.160 


4 


-0.6342 


4 


-0.6342 


3 


-0.6370 


3 


-0.6344 


0.180 


4 


-0.6671 


4 


-0.6670 


3 


-0.6678 


3 


-0.6672 


0.200 


3 


-0.6978 


3 


-0.6983 


3 


-0.6949 


3 


-0.6956 


0.250 


3 


-0.7682 


3 


-0.7689 


3 


-0.7475 


3 


-0.7493 


0.300 


3 


-0.8242 


3 


-0.8245 


3 


-0.7822 


3 


-0.7833 
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Fig. 3 - (Color online) Doping dependence of: (a,b) the magnetic incommensurability e, and (c,d) 
stripe filling v for the VSC (left) as well as VBC (right) stripe ground state deduced from figs. Qand 
121 for = (circles), t' = — 0.15t (squares), and t' — —0.3t (triangles). Solid line in (a,b) shows the 
idealized experimental behavior of e in LSCO [8]. 

Fig. 4 - (Color online) Doping dependence of the magnetic incommensurability e as obtained at 
U = 12t for: (a) (1,0) and (b) (1,1) spiral phase. Meaning of symbols as in fig. 01 



the system to maintain fixed filling v, but its actual value strongly depends on t' . This, in 
turn, results in almost doping independent chemical potential explaining the experimentally 
observed pinning of ^ [35]. In contrast, in the overdoped region with a lock-in effect of e, the 
size of the AF domains saturates, doped holes penetrate into the AF domains, and /i varies 
fast with doping. Unfortunately, for t' = — O.lSt, the established shift of exceeds the exper- 
imental value by a factor close to 2. Therefore, the present effective model can only explain 
qualitative trends and one needs to carry out calculations within more realistic multiband 
models including oxygen orbitals in order to reproduce quantitatively the experimental data. 

Finite t' < also modifies the spiral order — here one finds that the deviation of the 
spiral wavevector Q from the AF one strongly accelerates, as shown in fig. 01 It is worth 
noticing that the increase of the incommensurability with increasing \t'\ is consistent both 
with quantum Monte Carlo [36] and with DMFT [37] results. 

Summarizing, we have performed systematic studies of incommensurate phases by consid- 
ering variable size of AF domains in stripe phases, as well as the spiral order. Our findings 
obtained in the Hubbard model within the SB approach for t' = — O.lSt agree qualitatively in 
the low doping regime x < 1/8 with the experimental data for the cuprates and reveal a strong 
infiuence of the next-nearest neighbor hopping t' on the optimal filling of DWs. Simultane- 
ously, finite t' < promotes the onset of the diagonal spiral phases and increases the optimal 
spiral pitch e. We therefore conclude that a large value of \t'\ might be the reason why the 
static charge order has been detected in YBCO only in the highly underdoped regime [38]. 
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